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ABSTRACT 


This work utilized the computational techniques initially developed 
under NASA Grant No. NGR-47-005-085 and refined under NASA Contract No. 
NAS I — 1 02 1 0 , Task Order No. 2, to carry out further analytical studies 
of the transient and steady-state combustion processes in a hybrid rocket 
system. The particular system chosen consists of a gaseous oxidizer 
flowing within a tube of solid fuel, resulting in a heterogeneous 
combustion. Finite rate chemical kinetics with appropriate reaction 
mechanisms were incorporated in the model. A temperature dependent 
Arrhenius type fuel surface regression rate equation was chosen; for the 
current study. The governing mathematical equations employed for the 
reacting gas phase and for the solid phase are the general, two- 
dimensional, time-dependent conservation equations in a cylindrical 
coordinate system. Keeping the simplifying assumptions to a minimum, 
these basic equations were programmed for numerical computation, using 
two implicit finite-difference schemes, the Lax-Wendroff scheme for 
the gas phase, and, the Crank-Nicolson scheme for the solid phase. 
Appropriate initial and boundary conditions were represented, including 
heat and mass conservation at the gas-solid interface. A favorable 
downstream pressure gradient was imposed in order to obtain numerical 
solutions. In this study, Plexiglas (PMM) 1 and oxygen (C^) were 
selected as the fuel and the oxidizer, respectively. A single overall 
irreversible chemical reaction, such as fuel + oxidizer product, was 
used for the gas-phase chemical reaction to economize computing time. 
Constant property models were employed, using equivalent turbulent 
values for the transport properties. The initial and boundary conditions 
were chosen to simulate the conditions at the start of the ignition 
process when a fully-developed flow of preheated oxygen reaches the 
inlet of the combust ion chamber ; - Numerical calculations were performed 
for six cases with various values of flow rate and chamber pressure. 
Analyses were made primarily of the steady-state regression rate as a 
function of the axial position. The results were discussed. Conclusions 
have been drawn and recommendations have been made. 
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I. INTRODUCTION 


The work reported herein represents an extension of the work per- 
formed under NASA Contract No. NAS I — 1 02 1 0 , Task Order No. 2 [J], fol lowi 
initial work under NASA Grant No. NGR-47-005-085 [2j. The numerical 
computation was performed on the NASA Langley's CDC 6400/6600 computing 
faci I ity . 

The analytical models have been developed with a view to simulating 
realistic combustion processes in nonmetal ized hybrid rocket systems. 

In order to check out the numerical computational techniques used, 
previous calculations had been performed for hybrid rocket systems with 
axial lengths of 0.36 radius. The purpose of this investigation was to 
perform further tests on the analytical models, using a more realistic 
axial length of 3.7 radii. The temperature dependent Arrhenius type 
fuel regression rate model was chosen for the study. 

Again, the PI exig las-oxygen (PMM-O 2 ) system was selected for the 
current study. This allows the order of magnitude for the rate constant 
of the PMM.NQ 2 reaction to be estimated from the previous work in 
reference p], Parametric variations of flow rate and chamber pressure 
were made in an attempt to correlate the published experimental steady- 
state regression rate data. 

In the interior-bal I istic design of h^B.m'iidi rocket systems it is 
desirable to be able to predict the rate of fuel consumption as a 
function of axial position and of rate of oxidizer injection. The avail 
able steady-state experimental data were tailored towards this end. The 
steady-state regression rate as a function of axial position at apcnown 
oxidizer mass injection rate is, therefore, of prime interest in the 
current study. 

Numerical calculations were- performed for six-cases, and, profiles 

of the steady-state regression rate as a function of the axial position 
were obtained. Instead of an increase in the steady-state regression 
rate, there was a decrease when the pressure was increased at constant 
oxygen mass injection rate. Steady-state profiles of flame temperature 
and mass fraction of the product were analyzed. This adverse effect 



is believed +o be caused by a layer of cold unburnt mixture of fuel and 
oxidizer in the immediate vicinity of the fuel surface acting as a heat 
shield for the fuel surface and thus blocking off the increased heat 
release from the combustible gas at higher pressure. Recommendations 
have been made for corrections to vary the chosen inlet conditions and 
the values for the reaction criteria to eliminate this effect. 

It is shown that the flame position in the combustion chamber can 
be predicted from the calculated results, and that, this type of infor- 
mation can possibly be used for future experimental correlations. 
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II. ANALYTICAL MODELS 

The system under study, as shown In Fig. I, consists of an oxidizer, 
such as oxygen, flowing through a tubular solid material where chemical 
reactions take place in the gaseous phase. The fundamental approach has 
been followed in formulating the analytical models. The mathematical 
equations derived have been expressed in dimensionless form. Finite 
rate chemical kinetics and appropriate reaction mechanisms have been 
incorporated in modelling the gas phase chemical reaction. A temperature 
dependent Arrhenius type fuel regression rate model has been selected 
for the current study. 


r 
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A. Non-Dimensional Equations; 


The mathematical formulation employs the general time-dependent 
two-dimensional conservation equations in cylindrical coordinates as 
the governing equations for the gas phase, the solid phase, and, the 
gas-solid interface £2 J. The governing equations: in the gas phase are 
as fol lows : 


Species Continuity Equations : 

A A 

- Dy I 3 3y.| PD. 9y . 9 ( ^ 9y ; 

p' —r. — = C - CpD. /s ) + ( ^ + —s- (pD-. — * — )] + DaW. , 

DDt Re Sc 9r 1 9r r 9r 9z 9z 1 

(I ) 

i = I , 2, . . . , N 


Momentum Equations: 
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Energy Equation: 
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( 4 ) 


♦ v A A 

t ®-l DgCS h.w.) 

.. ii 
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Equation of State: 


P = iCZ (ViV^T 


N 


Ec. 


i = l 


The governing equation in the solid phase is as follows: 
Transient Heat Conduction Equation : 


(5) 


P 


C 


s s 



so 


F V 

o o 





V X 9T 9 . 9T ■> 

(-^-<^J0-+ ( rr4\ l -Af] 
r 9r 3z s 3z 


(6) 


The governing equations derived for the interface are as follows: 
Mass Balance: 


D 3y . 
o , 1 1 


r v > ~ 

rovo 3r w 
o o 


p . D. 
w i w 


- p v y. 
w w' i w 


p so1 , r o* • 

( — } ( - )p s r y isw’ 
p o V o 


i = I ,22, . .., N 

Energy Bl 
Energy Balance: 


(7) 


T a.) X a * 3T . . r v ■ i * aT 

SO SO Q ^ 1 o «v 0 1 

‘r 2 ’ (x s 7* • v, <— > B <(.,», >] - 

T o X o; 8r " “» 1 3r w 


T . X r r 

+ P s / V 2 -) (h .y. )] 

SW T o X o °ioor i sgi y ,sw J 


Bo - - 3y i 

p w #) ci (h.o. ^J.) 


( 8 ) 


ct j / l 
o 


3r w 


By, definition, the sum of the mass fraction of chemical species is 
equal to unity. Hence, there are only -(N-l) independent species 
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continuity equations. However, in order to e I imi nate roundoff error 
in the process of subtraction by computers, N of the species continuity 
equations were solved. 

The complete set of equations should include the continuity equation, 
however, difficulties with the numerical solutions were experienced when 
it was included. As an alternative, the set of equations shown above 
was solved with an initially assumed pressure gradient. Numerical 
iteration of the pressure gradient should then be performed to obtain 
the exact solution. 

B. Numerical Methods Used_ 

~ — ■ ■ 

The appropriate finite-difference schemes of second order accuracy 
have been selected to solve these dimensionless equations numerically. 

An implicity Lax-Wendroff scheme for the gas phase equations, an implicit 
Crank-Ni co I son scheme for the transient heat conduction equation, and a 
time-independent finite-difference scheme for the boundary equations have 
been employed in deriving the difference equations for computer pro- 
gramming. The difference analogs of these fiit'nite-d i fference schemes are 
presented in reference [3]. 

Since the numerical schemes employed are implicit, an iterative 
procedure had to be used. In -addition, a telaxaf ion/procedure was 

-- S'. . 

employed for the purpose of convergence, and, the A itken’s -6.2/ me +hod [43 
was used as an accelerating convergence procedure during the iterations 
at each time level to speed up the calculations. 

In the current study, the combustion process is dominating. An 
imposed constant favarable downstream pressure gradient was found 
sufficient. Hence, no numerical iteration of the pressure gradient has 
been performed. 

C. Reaction Kinetics 

Since there is a lack of understanding of the kinetic mechanisms 
for most of the f ue I -oxi di zer combinations of interest, a one-step stoi- 
chiometric reaction is ! assumed as follows: 
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( 9 ) 


Vf I _ + v-. I~ v n I 0 
F F Ox Ox P P 


The reaction mechanisms are modeled by two criteria - a critical temp- 
erature, T , and, the upper and lower limits of flammability (UFL,LFL). 
Mathematically, these two criteria are as follows: 


T > T 


( 10 ) 


UFL > LFL 

Vox 


(II) 


When the above two criteri a|are met locally, the reaction is assumed to 
occur. The rate of production of species for this overall reaction is 
calculated as follows: . 


: k i; M. 

w. = SLL 

I W 


(V i y . ' ),Ej p3a Ej n 3o 

1 liQxj 


(12) 


where a and b, are the react ion . orders for the fuel and ox'idfe. 

where a and b are the reaction orders for the fuel and oxidizer respectively. 

D. Regression Rate Equation ■ 

It has been demonstrated experimentally that the fuel surface re- 
gression rate in hybrid rocket systems varies in the axial direction. 

The regression rate model should allow r to be calculated from the input 
parameters and be used as a combustion model to predict trends in r 
due to operational variables. Hence, the current model assumes that . 
the fuel surf-ace ..chemical'' reaction; i s .a - f'unct Lon.-of -the local, temper-' 
at.ureeatjatheogas'TSO I Tdf JmtserfaceyuTh *tsuch •: thatcthetf uetjsurface 

regression rate can be expressed in the Arrhenius form as follows: 

.'regress 


r = A s exp 


EE \ 

tss 

R T“ 


( 13 ) 


where A g and E g are empirical constants. 
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III. NUMERICAL CALCULATIONS 

The hybrid rocket motor selected for the current analytical study 
has the following dimensions: r = 2 .54 cm ( I .0 i n ) , r = 2.56 cm (1.007 in), 

and, L = 9.408 cm (3.70 in). 7 For the current study, six sets of 
ca I cu I at ionswwere performed with pre-selected values of flow rate, and 
chamber pressure. 

A. Selection of Fuel and Oxidizer 

Plexiglas (PMM) and oxygen (C^) have been selected as the fuel and 
the oxidizer for the current study ,t> because our previous work £03, ;a I lows 
us to pick the correct magnitude for the rate constant in mode I I i ng the 
PMM-O 2 reaction through parametric variation of the rate constant. In 
addition, the PMM-0 ^ combination is a relatively simple system which 
can readily be simulated in an experimental program. 

B. Initial and Boundary Conditions 

The initial and boundary conditions used were as follows: 

(a) A fully-developed flow of pre-heated oxygen (1620° R) reaching 
the inlet of the combustion chamber, such that the axial component of 
the velocity distribution at the inlet is based on the (l/7)-power law 
for turbulent flow, i.e., 

" "1/7 

v z = (I - r) 1 ; (14) 

(b) The assumed initial temperature profile represents the conditions 

at the start of the ignition process; 1 

(c) The temperature at the exterior surface of the solid fuel is 
constant ; 

(d) The axial component of velocity is zero at the gas-solid inter- 
face and the radial component of velocity is zero at the center line. 
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C. Reaction Kinetics 


Since there is a lack of understanding of the kinetic mechanism 
for the PMM-O^ reaction, a one-step sto i ch iometr ic reaction is assumed 
as foil ows : 


C 5 H 8 0 2 + 60 2 ^ 5C0 2 + 4H 2 0 (15) 

In order to save computer time, the reaction products C0 2 and H 2 0 are 
treated as one single species such that the reaction can be represented 
as fol lows: 

C 5 H 8 0 2 + 60 2 & 9 Product (16) 

where the product is a mixture containing 5 parts of gaseous carbon 

V 

dioxide and 4 parts of water vapor. 

Treating the overall reaction of PMM and 0^ in the same manner as 
in reference PH, an order-of-magn i tude value for the rate constant was 
estimated for the current study. This selected value was _ ,~'j 
k = 3.9 x 10^ cir/Vg-mo I e 2 /sec (10*^ ft^/ 1 b m ~mo le^/sec) . 

D. Property Models 

The models for the transport properties, the fluid properties and 
the properties for the solid fuel were assumed to be constant in the 
same manner as reported in reference [2]. The property values were 
estimated if they were not readily available. The selected thermophys ica 
properties for PMM, 0^ and the product were obtained either directly or 
by estimations from references PH through [9] and are tabulated in 
Table I. Following the approach outlined in reference [2D> equivalent 
turbulent values were used for the transport properties. These selected 

7 

values are y = 0.00893 g/cm/sec (0.0006 Ib^/sec/ft), D. = 8.23 cm /sec 
(0.00886 ft 2 /sec ) , and X = 0.0031 ca I /sec/cm/°K (0.75 Btu/hr/ft/°.R) . 



TABLE I 

Selected Thermophysical Properties for Plexiglas (PMM), 
Oxygen and the Product* 



PMM 

°2 

Product*..; 

g/g-mole, or, 

v 

lb /lb -mo 1 e 
m m 

100.0 C5U 

32.0 [5J 

32.4 [5] 

cal/g/°k, or, 
Cp, 

Btu/lb /°R 
m 

0.35 [6] 

0.25 [7] 

0.20 [7] 

ca 1 /g 
h°, 

Btu/ 1 b 

m 

- 928 
-1670 [8] 

0 

62260 1 ", 7 
-4680.0 [5] 

°K 

Boiling poi nt, 

°R 

374 

673.2 C8] 

90. 19 
162.34 [5] 

373.3 
672.0 [5] 

, t ca 1 /g 

si, Btu/lb 
’ m 

130 

234 [8] 

- 

- 

, tv caj/nA 
btt ca l//g 

■tig , Btu/ib m 

863 

86 | 
155 [8] 


r.r 3 3.4 

240.0 E5] 

g/cm 3 

p Q > 

lb /ft 3 
m 

1.683 
150 [6] 

- 

- 

cal /sec/cm/°K 

X , 

S Btu/hr/ft/°R 

0.0012 
0.289 [6] 


- 

ca 1 /g/°K, or, 
C , 

S Btu/lb /°R 
m 

0.31 [9] 

. . 



*A mixture containing 5 parts of gaseous carbon dioxide and 4 parts of 
. water vapor. ^ 

'heat of fusion heat of vaporization 
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Figure 2. Grid Layout 


E. Imposed Pressure Gradient 


In order +o obtain numerical solution, a constant favorable down- 
stream pressure gradient equal to that used by Jones and Issaacson [10] 
in their analysis was imposed as in reference [2]]. This imposed 
pressure gradient is 

- ^ = 2.576 dyne/crn /cm (0/1.64 -psf/f« . (17) 

F. Grid Layout 

The grid layout selected for the numerical calculations is shown 
is Fig. 2. There were 15 stations along the axis in the ;z'< di recti on. 

The ratio of the grid size for the three uniformly spaced regions in 
the axial direction was 1:2:4. "In the r direction, there were 20 
sitiatliisnsatnpitrheegasdpliiase and 10 stations in the solid phase. Finer 
grid size (7:1 for the gas phase and I ?5 for the solid phase) was 
employed near the gas-solid interface. 

G. Empirical Constants for the Regression Rate Equation 

A trial and error method as discussed in reference ,C! 3 was used 

to determine the values for the empirical constants A g and E g such 

that r = 0.254 mm/sec (0.01 in/sec) at T = 644.4° K (1160° R) , the 

w 

initially assumed temperature value at the fuel surrfjaee. Tlhese-; __ 
numerically determined values are 

A g = 0.645 mm/sec (0.002116 ft/sec) (18) 

E = 1.0 Kcal/g-mole (1799 Btu/lb -mole). (19) 

s ■ m 

H7 Parametric Values Used ' " * ~ 

Numerical calculations were performed for six cases with various 
values of flow rate and chamber pressure. A compendium of the parametric 
values used for the numerical calculations of the six cases is shown 
i n Tab I e II. 


12 



Table II 


Parametric Values Used for the Numerical Calculations 


1 

Case 

V 

vo^ 

o 

G 

c 

1 

) 

P 

i 

m/sec 

ft/sec 

g/cm^/sec 

2 

lb /in /sec 
m 

atm 

1 

32.92 

1080 

18.98 

0.27 

1 .67 

2 

36.58 

120 

6.32 

0.09 

5.0 

3 

7 

360 

18.98 

0.27 

5.0 

4 

329,2 

1080 

56.94 

0.81 

5.0 

5 

36.58 

120 

18.98 

0.27 

15.0 

6 

|09.7 

360 

56.94 

0.81 

15.0 








IV. RESULTS AND DISCUSSIONS 


In performing the numerical calculations, both, the transient 
and the steady-state results have been obtained, the transient 
solutions provide the guideline in determining when steady-state 
conditions have been reached. In general, the total time necessary 
for the transient solution to reach steady-state is approximately 
twice the time required for the oxidizer at the inlet to reach the 
axial station being considered. 

Since no transient experimental data on this type of system 
is available, analyses were made primarily on the steady-state results. 
Hence, only the steady-state solutions are included to illustrate the 
combustion induced variations in flame temperature, species mass 
fraction, and regression rate profile as a function of input parameters. 

In all the six cases calculated, the initially assumed profiles 
for the dependent variables, v^, v z> T, T s , p, and, Yq at any 

axial station are shown in Figures 3a-g, respectively. I n it? ally, 

the mass fraction of thep'prdductv/wasaassumedttobbezzenoeeverywhere. 

Holding the pressure constant at 5 atmospheres, the steady-state 
regression rate for cases ^2.?/ 3 and 4 (see Table II) are plotted in 
Figure 4. It can be observed that an increase in flow velocity at 
constant pressure, or, in other words, an increase i n the mass flow 
rate of oxygen at the inlet tends to increase the steady-state regression 
rate. Similar trends were observed at a constant pressure of 15 atmo- 
spheres, i.e., cases 5 and 6. 

The steady-state regression rate was found to be independent of 
pressure at flow velocity values of 109.7 m/sec (360 ft/sec) and 
329.2 m/s ec (1080 ft/sec). This is illustrated by plotting the steady- 
state regression rate of gases 3 and 6 at constant flow velocity of 
109.7 m/sec (360 ft/sec) as shown in Figure 5. However, at lower flow 
velocity of 36.58 m/sec (120 ft/sec), an increase in pressure tends- 
to increase the steady-state regression rate. This is demonstrated in 
Figure 6. 
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With the oxygen mass injection rate held at a constant value of 
2 2 

18.98 g/cm /sec (0.27 I b / i n /sec), the steady-state regression rate 

profiles of cases I, 3 and 5 are plotted in Figure 7. It is observed 

that an increase in pressure at constant oxygen mass injection rate 

tends to decrease the steady-state regression rate. The same behavior 

was observed for cases 4 and 6 when the oxygen mass injection rate was 

2 2 

held at a constant value of 56.94 g/cm /sec (0.81 lb /in /sec). How- 

m 

ever, Marxman, Woolridge and Muzzy p I H, in their experimental study 
of PMM-C^ systems, observed an increase in the steady-state regression 
rate at higher pressure. In addition, this same effect was reported 

S 

by Smoot and Price p2Q| in their experimental study of nonmetal ized 
hybrid fuel systems using three different binder compounds (ButyL/ 
rubber, PBAA and Po I yruethane ) with oxidizers varied from 100$ fluorine 
to 100$ oxygen. Hence, the current analytical results do not confirm 
their experimental observations.. The scattered experimental data 
shown in Reference P I H for PMM-C^ systems at atmospheric pressure are 
shown in the cross-hatched region in Figure 8. For comparison, the 
steady-state regression rate values at axial station 14 for the six 
cases ca I cu I atedjhare plotted in Figure 8. 

Steady-state profiles, of flame temperature and mass fraction of 
the product £f'o r the six cases calculated were reviewed. It was observed 
that an increase in pressure at constant oxidizer mass injection rate 
did increase the flame temperature considerably as shown in Figures 
9a-b, and the mass fraction of the product as shown in Figures lOa-b. 
However, from these figures, it was also observed that the gas-phase 
combustion did not sustain itself at the radial station next to the fuel 
surface. This layer of cold unburnt gas mixture decoupled the heating 
effect on the fuel surface from the combustion zone. As a result, the 
calculated steady-state regress ion- rates d i d not display the correct 
kinetically controlled effect, tie. , an increase in the steady-state 
regression rate with an increase in$pressure. 

With the proper choice of'/the profiles of the mass fraction of fuel 
and oxidizer at the inlet such that it is less fuel rich at radial station 



next to the fuel surface, and, with values for the reaction criteria such 
that combustion can sustain itself at the radial station next to the fuel 
surface, it is believed that the correct pressure effect at constant 
oxidizer mass injection rate should result. 

The current analytical method also provides the capabi I ity of 
predicting the position of the flame in the combustion chamber. To 
demonstrate this point, the peak of the temperature profile was taken 
as the flame position, and, the steady-state flame position in the 
combustion chamber for case 3 was predicted as shown in Figure II. 

Using Schlieren and shadowgraph techniques, Muzzy [J3H was able to 

-v, 

predict experimentally the steady-state flame position in a slab burner. 
The analytical prediction is in qualitative agreement with Muzzy ' s 
experimental results. Hence, it is believed that this type of technique 
can possibly be used for future experimental corre I at ions . 
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V. GONCLUS IONS AND RECOMMENDATIONS 

The computer program developed under NASA Contract No. NAS I -1 02 10, 

Task Order No. 2 following initial work under NASA Grant No. NGR-47-005-085 
for the study of transient combustion processes in hybrid rocket systems 
has been modified to allow a study of systems with a more realistic 
axial length (3.7 r ). Current emphasis has been placed on correlating 
the calculated steady-state regression rate results with published 
experimental data. The following has been achieved: 

(a) The computer program can now handled a more realistic axial 
length (3.7 r ) with reasonable usage of computing time. 

(b) The steady-state regression rate at constant pressure increases 
with increased oxidizer mass injection rate. 

(c) The computed steady-state regression rate at constant oxidizer 
mass injection rate decreasedwwtthTiiincre'ase'dppeessuBewwIii&h i i;s in .const rast 
with reported experimental observations. More work is required to 
investigate these results, which appear to be incorrect. 

(d) The steady-state results can possibly be used to predict the 
flame position in the combustion chamber. 

A list of recommendations is as follows: 

(a) The inlet mass fraction of fuel and oxidizer at radial 
station next to the fuel surface and values for the reaction criteria 
should be cyaried such that combustion can sustain itself at radial 
station next to the fuel surface further downstream. 

(b) Concurrent analytical and experimental work should be performed 
to confirm the validity of the current analytical models. 
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FIGURE ( 9a ) 
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